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Abstract. The asymmetric simple exclusion process with random-force disorder is 
studied within the mean field approximation. The stationary current through a domain 
with reversed bias is analyzed and the results are found to be in accordance with earlier 
intuitive assumptions. On the grounds of these results, a phenomenological random 
barrier model is applied in order to describe quantitatively the coarsening phenomena. 
Predictions of the theory are compared with numerical results obtained by integrating 
the mean field evolution equations. 
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1. Introduction 

Transport processes in nature, like intracellular transport which is reahzed by active 
motor proteins [1] are often modeled by simple exclusion processes, in which particles 
residing on the sites of a lattice hop stochastically to neighboring sites provided the 
target site is empty [21 EJ IH El E]. For this paradigmatic model of driven interacting 
particle systems many exact results are available [7] . As most of the real systems are not 
ideally translationally invariant, for instance the filaments on which molecular motors 
move are heterogeneous, a challenging problem is the study of spatially inhomogeneous 
versions of exclusion processes P|9l[Tni[IIl[l2l[l3l[Tlll5l[T6l[I71[T8]for which the bulk of 
results is obtained by phenomenological methods based on the statistics of extremes, by 
mean field approximation and by Monte Carlo simulations. Most works concern the one- 
dimensional totally asymmetric process where particles can hop only in one direction 
with site dependent quenched random rates. In such systems clusters of consecutive 
bonds with low hop rate act as bottlenecks and the stationary state is segregated, i.e. 
consists of macroscopic regions of low and high density [10]. When the system is started 
from a state with homogeneous density it undergoes a coarsening process in which the 
typical size of low and high density segments is growing in time [TOl US]. A similar 
coarsening phenomenon occurs in the partially asymmetric simple exclusion process 
with random- force disorder, where the direction of the local bias is random [ini [131 US] ■ 
In this case, clusters of bonds with reversed bias compared to the global one limit the 
current and since their extension is unbounded in an infinite system, they result in that, 
parallel with the coarsening of the length scale, the local currents tend to zero in the 
long time limit t — t- cxd [15]. In the driven phase of this model, a phenomenological 
random trap description was developed which relates the coarsening exponents to the 
dynamical exponent of random walk in random environment and the predictions of this 
theory has been found to be in agreement with Monte Carlo simulations [15]. 

In this paper we shall investigate this model within a mean-field approximation, 
which, to our knowledge, has not been applied to the disordered partially asymmetric 
model yet. Calculating the steady state current through a segment with a reversed bias 
we shall argue that, in the driven phase, extreme value statistics of barrier heights leads 
to the same dynamical exponents in the less complex mean field model as those of the 
original one. As opposed to earlier works applying mean field approximation, here we 
focus on the dynamical behavior rather than the steady state. The phenomenological 
predictions will be checked by numerically integrating the dynamical mean field 
equations. 

The rest of the paper is organized as follows. In Sec. |2] the model is defined in 
details. In Sec. |3l the elements of the phenomenological theory of the steady state are 
surveyed. Sec. Hlis devoted to the analysis of the current through a single barrier within 
the mean field approximation with different boundary conditions. The phenomenological 
theory of the dynamics is reviewed and applied to the model in Sec. |5]and the predictions 
are compared with numerical results in Sec. [61 Finally, the results are discussed in Sec. [7] 
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and some calculations for the dynamics of the pure model are presented in the Appendix. 
2. The model 

The disordered partially asymmetric simple exclusion process is defined as follows. An 
infinite one-dimensional lattice is given the sites of which are either empty or occupied 
by a particle. On this state space a Markov process is considered in which particles hop 
independently to an adjacent site provided that site is empty. The hop rate from site 
i to site 2 + 1 — 1) is denoted by pi (qi-i) and the pairs of rates {pi, qi+i) are i.i.d. 
positive random variables. Furthermore, we require that < Prob(pj < qi-i) < 1. In 
words, the local force Fi = ln{pi/qi) acting on particles can be both positive or negative 
with finite probabilities. 

In the mean field approximation, the pair correlations of the occupation number 
Hi = 0,1 are neglected meaning that expected values of products of occupation numbers 
{uiUi^i) are replaced by (nj)(nj_|_i). Then the evolution equation for the local density 
Piit) = {ni{t)) reads as 



Defining the model on a finite ring of sites L rather than on the integers it has a 
steady state where the local currents Jj are all equal. This stationary current is sample- 
dependent i.e. depends on the set of random hop rates {pi,qi}. The typical stationary 
current in the ensemble of samples of size L tends to zero in the limit L — t- oo due to 
the occurrence of larger and larger domains with reversed local force that control the 
current [TOl |T5]. In the infinite system, the local densities do not converge in the limit 
t — )■ oo therefore there exists no stationary state. Nevertheless, when the system is 
started e.g. from a homogeneous state, the local currents Ji{t) which are non-zero for 
finite times all tend to zero in long time limit [15]. We shall consider the dynamics of 
this non-stationary process and are mainly interested in the dependence of the typical 
current 



Before analyzing the disordered model, we discuss the evolution of the typical 
current in the pure model where Pi = p, qi = q for all i. As it is shown in the Appendix, 
the typical deviation of the current from the stationary one (Joo) decays algebraically 
with the time. The decay exponent depends on the symmetries of the model. If p = q 
(symmetric simple exclusion process) the typical current decays as 






on time. 



Jtyp{t) ~ t 



(4) 
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U p ^ q (asymmetric simple exclusion process) then, for densities different from 1/2, 
the typical current decays as 

(J- Joo)typ(t) (5) 

while at half- filling p = 1 /2 we have 

{J - Joohypit) ^ t-^/\ (6) 

These latter results follow essentially from the time-dependence of the typical deviation 
of the local density from the stationary value calculated by Burgers [19] but for the sake 
of self-containedness a short heuristic derivation is given in the Appendix. 



3. Phenomenological random barrier theory 

For exclusion processes with random-force disorder a phenomenological theory exists 
by which many steady state and non-stationary properties are successfully described in 
accordance with results of Monte Carlo simulations [TOl [15]. The basic idea is that the 
random environment (i.e. the series of jump rates) contains localized trapping regions or 
barriers in which the local force is reversed compared to the global one and such regions 
therefore can maintain a very low current. These barriers can be defined quantitatively 
in terms of the potential Ui which is defined by 

AUi = U,+i -Ui = -F, = \n{q,/pi). (7) 

An interval from site a to site b is said to be an ascending interval if and Ua < Ui < Ub 
for a < i < b. The ascending interval [a, b] is a barrier if there does not exist a longer 
ascending interval which contains [a,b]. If the average force is non-zero, i.e. AU ^ 0, 
where the overbar denotes averaging over the distribution of hop rates, the system is in 
the driven phase and the size of the barriers has an exponentially decaying distribution 
and the number of barriers in a finite system is proportional to the size of the system. 
Each barrier has a maximal carrying capacity and the smallest one among these values 
determines the stationary current of the finite system. The key question in this theory 
is how the maximal carrying capacity varies with the parameters of the barrier, which 
depends on the particular model. In case of the partially asymmetric simple exclusion 
process this has been obtained by the following phenomenological arguments. The 
steady state of a homogeneous, open system with reversed bias (g > p) where particles 
enter at site 1 with rate a and are removed at site L with rate /3 is exactly known |20j . 
The density profile contains an anti-shock in the middle of the system, which separates 
a high density phase on its left hand side where the density is close to one from a low 
density phase on its right hand side where the density is close to zero. In case of an 
inhomogeneous barrier the exact steady state is no longer available but the profile is 
qualitatively similar to that of the pure case. In the steady state, the anti-shock must 
be located where the potential (measured from the bottom of the barrier) is half of the 
total height of the potential, since the current of a single particle in the low density 
phase must be equal to the current of a single hole in the high density phase. Since the 
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distribution of heights of barriers can be calculated, the distribution of the current in 
finite systems is obtained by applying the statistics of extremes [15], [16]. 

We will apply this phenomenological theory to the mean field model defined above. 
First we calculate the mean field current through a random barrier and shall see that it 
is determined practically by the height of the barrier as it has been intuitively assumed 
for the original stochastic model. 



4. Mean field current over random barriers 



4-1- Open boundaries 

Let us consider an open random barrier with sites and with entrance and exit rates a 
and /3, respectively. In the steady state, the local densities in the bulk obey the relations 

PiPi{l - pi+i) - Qipi+lil - Pi) = J, (8) 
where the current J is to be determined. Introducing the variables = pi/ (1 — pi) and 
=Pi/Qi Eq. ([8]) takes the form 

Vi+i = riVi - Jql^{l + yi){l + Vi+i). (9) 

This is a non-linear recursion equation for the densities. Let us choose a site where the 
potential measured from the left end of the system is roughly half of the total height of 
the potential barrier and denote this site by 0. As we shall see later this site is in an 
anti-shock region where the density is close to 1/2. Thus yo ~ 1. Denoting the term in 
Eq ([9]) which is responsible for non-linearity by 

UJ^ = qr\l + yi){l + y,+i) (10) 

and regarding it as if it was a constant, the recursion can be formally carried out 
starting from site to the right, i.e. toward the low density phase till the rightmost site 
L, yielding: 

/L-1 \ r L-1 j 

yL=[\{rA yo-J Y^^^H^^^ 

\j=0 J I j=0 i=0 

The current J is simply related to the density at this site as follows: 

J = f3pL = /3t^ = /32/L + 0{yl). (12) 

1 + 2/L 

Here, we have used that, as we shall see a posteriori, the density pi, as well as yi decay 
exponentially with the site index, and they are thus very small for large L. Eliminating 
yi from the latter two equations, we obtain the following formal expression for the 
current: 

L-1 L-1 

AL = J2^,r-'l[r,. (14) 

j=0 i=j 



with 
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The variable A^^ is a function of the densities pi but, as we shall show below, it is bounded 
by a random variable which is finite [0(1)] in typical barriers. Since the second term in 
the brackets on the r.h.s. of Eq. ( |TT|) is negative, the inequality 

i-i 

Vj < yoYlri = yoe~^' (15) 

1=0 

obviously holds for all j > 0. Here and in the following, the potential at site is set to 
zero, i.e. Uq = 0. Using these inequalities, we can write 

L-l 

0<Al< e-^^ 5^g7^(e^^+i + yor;^){yoe~^'+' + 1) = 

j=0 

As can be seen, only those sites give an 0(1) contribution to this sum at which the 
magnitude of the potential is close to Ul, i.e. either Uj ~ Ul or Uj ~ —Ul. The main 
contribution comes from the sites in the end region at which the potential is close to 
Ul- Since the random potential is, in general, not necessarily monotonic there may 
be also sites far from the end with \Uj\ ~ Ul- Nevertheless, the barriers have a finite 
(non- vanishing) average slope Ul/L in the limit L — )■ oo, therefore the number of such 
sites, as well as the random variable A^, is expected to have an L-independent limit 
distribution. If the potential does not turn down to the vicinity of —Ul-, which is the 
typical situation, we can obtain an accurate estimate of A^ as follows. In this case 
the terms yi and appearing in Eq. (JH]) through Uj can be neglected according to 
inequality f|T5|) . This results in the following expression 

L-l L-l 
i=o i=j 

= Ql-i + C~2^L-i + (lL^3^L~irL~2 + --- + qo^rL-irL-2 ---n, (17) 

which is thus an accurate estimate of A^ for large L in the case of barriers for which 
the potential is well separated from —Ul, i.e. —Ul ^ Ui- Moreover, this sum starting 
with the term j = L — 1 as written above is rapidly converging if the potential is well 
separated also from Ul apart from the region close to the end of the system. This is 
the case for barriers with monotonic potential. In general samples or for finite L, A° is 
a lower bound on Al- 

Notice that neglecting the terms y^ and in Eq. (|9]) results in a linear recursion 
which describes independent random walkers with density yi at site i- As a consequence, 
the sum A^^ can be related to properties of random walks, as follows. Let us consider 
a finite lattice with sites 0, 1, . . . , L, L + 1 and the same series of hop rates as given for 
the exclusion process except that po is set to zero, furthermore Ql = and pl = P- 
That means, sites and L + 1 are absorbing. Starting at site 1, the probability that 
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the walker is absorbed at site L + 1 when t — oo is called persistence probability and 
is given by [2T| : 

L-l i 



Ppcrs (i) 



i=i j=i 



This can be recast as PpersiL) = Pq ^e^^^(A° + (3^^)^^ which leads to that, whenever 
the replacement of A^, by A° is justified, the current is asymptotically proportional to 
the persistence probability of the corresponding random walk: 

J° = yoPoPpcTsiL). (19) 

The expression of the current in Eq. (I13p is still incomplete in the sense that it 
contains the variable yo at the chosen reference site in the anti-shock region. This can 
be, however, easily eliminated as follows. Introducing the variables Xi = y^^, one can 
write the recursion in Eq. ([8]) for decreasing indeces in the following form: 

Xi-^i = ri^iXi - Jqi\{l + Xi){l + Xi^i), (20) 

which has the same structure as Eq. ^ for the forward iteration. Performing the 
recursion from the same initial site as for the forward iteration to the entrance site 
indexed by — L', and using the relation between X-l' and the current: 

J = a{l- p^L') = ax^L' +0{x\,), (21) 

we obtain an expression for the current analogous to Eq. (fT3!) : 

J = ^0 _f'' + 0{x\,), (22) 

where 

-L'+l 

^-L' = e""-^' Yl &^ + ^^^(l + (23) 

j=0 

Here, A^^/ has the same properties as A^,, e.g. the linear contribution for barriers with 
Uj -C —U^L' is given by 

A° = qZ\,+qZ\,+ir^L'+qZu+2r-L'r-L'+i+- ■ ■+qZ\r-L'r-L'+i ■ ■ ■ r-2-(24) 

This can be again related to a persistence problem in a finite system with sites 
—L' — 1, — L', . . . , —1, and with hop rates p[ = qi, q[ = pi for i = —2, —3, . . . , —L' 
and q'_i = pZy^i = 0, p'_i = q'_j^i_i = a. Now, the walker starts at site —1 
and the probability that it ends up at site —L' — 1 can be written as Pp^j.^{L') = 
(p_i)-^e^-i'(A°^, + a-^)-^. This leads to J° = Xop_iPp^,^{L'). 

Obviously, the current in Eq. (12^ must be equal to that obtained by the forward 
iteration in Eq. f fT3|) . Multiplying the right hand sides of the two equations and 
introducing the total height of the barrier as Un = Y[j=-L' = Ul + \U^l'\, we obtain 
the following formal expression for the stationary current: 

J = [{a-' + A^L'){/3-' + A^)] e-^-/2 + 0{J^). (25) 
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Although Al and A_i/ in this expression are given in terms of the density profile {pi} 
which is not known exactly in a closed form, they are bounded by random variables 
which are typically 0(1). Moreover, if relations 

Ui > U^L' for i > 0, 

f/, <^Ul for i < (26) 

are satisfied then A^^ and A^^/ are accurately approximated by the linear contributions 
given in Eqs. ( ITTl) and ( 12^ for large barriers. In this case, the current can also be given 
in terms of persistence probabilities of random walks: 



J' = v/prip^y'ppers(i^)Ppers(^0- (27) 

In fact, it is easy to see that the only barriers for which the above approximations are 
invalid are those which have a bulk site with Ui ~ Ul and another one with Uj ~ — 
furthermore i < j. In all other cases the reference site to which the summations in 
Al and A^^i go, can be shifted such that sites with Ui ~ Ul {Ui ~ —U_l') are on the 
right (left) hand side of the reference point. 

In case of a homogeneous barrier with Pi = p and qi =^ the condition in Eq. fl26!) 
is obviously fulfilled and the current for large is given b} 

af3{q — py 



AN) 



r 



(N-l)/2_ .28) 



(a + q-p){P + q-p)_ 
For this asymptotically exact mean field current an approximate formula has been 
derived in Ref. |H] in the limit of weak asymmetry p ^ q. 



4.2. Barrier in an infinite system 

Regarding that the disordered model contains random barriers embedded in it, we will 
consider boundary conditions more appropriate for the above problem, namely the 
maximal current through a single barrier which is part of a large disordered system 
will be analyzed. 

Let us assume that the barrier is very far from other barriers, i.e. the potential is 
monotonically decreasing outside the barrier. Starting the iteration again from a site in 
the barrier where Uq ^ 1, we have 

l/n = e-^"(yo-^e^"7n), (29) 
where we have introduced the variables 

n-l 

7. = J]a;,e^^+-^^ (30) 

j=0 

Outside the barrier, the mass flows with a non- vanishing 0(1) velocity and taking into 
account that the stationary current is 0(e~^^) this implies that ?/„ = 0(e~^^) for n ^ L. 

I The exact current of the asymmetric simple exclusion process calculated in Ref. |20) differs from this 
mean field current by a factor of r^/^. 
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Thus, far from the barrier ?/„e'^" tends to zero and, in the hmit n — )■ oo, we obtain from 
Eq. m 

J ~ yo (31) 

Too 



for large L. A similar backward iteration yields J ~ Xq^-t— where 'j'_oo = + 
Xj){l + Xj-i)e^-L'~^^'^ . Multiplying the two expressions for the current yields finally 

J^{looi.^r"'e-^-'\ (32) 

The sum can be decomposed as 

L— 1 oo oo 

^^ = Y^ w.-e^^+i-^^ + a;,-e^^+^-^^ = + ^ a;,-e^^+i-^^ (33) 

i=0 j=L j=L 

In the limit of large barriers {L — )■ oo) we have coj ~ qj^ for j > L and thus 

7oo ^ Az. + + Pll^rl^ + PlWlliTl' + ... (34) 

The first term on the r.h.s. is the same one that appears in the current of an open 
barrier while the sum of the other terms converges since outside the barrier the potential 
decreases monotonically (rj > 1). In case condition ( 126|) is met such as for monotonic 
barriers, the term A/, in the asymptotic form in Eq. ( l34l) can be replaced by A^^. 
Furthermore, in this case 7oo is related to the persistence problem in a semi-infinite 
lattice with sites 0,1,2,... where the walker starts at site 1. The probability that, for 
t — )■ oo, the walker is not at site can be written as Ppers(oo) = \\mL^ooPpcvs{,L) = 
(Poe^-7i)-^ 

The sum 7^^^ can be written in a similar form: 7^^^ ~ ^-u + Pli'-i + 
pZif _2''"~L' -1 + P~L'-3^-L'-2'^-i'-i + • • • for the linear contribution we have the 
relation Pp^,^{oo) = limi/^ooPpers(^') = [P-ie~^-^' (7-oo)T^- 

Let us now consider a homogeneous barrier where Pi = p and qi = q ii —L' < i < L 
and Pi = q and qi = p otherwise. In this simple case the asymptotic forms of 700 and 
jLoo can be easily evaluated yielding the current for large A^: 

j(iV) ~ i_Z^(A'-i)/2. (35) 
4-3. A more general model 

As can be seen in the form of the current in Eq. ([2]) , the factors of the form 1 — pi ensure 
that the local density at any site, provided it was initially below 1, cannot exceed this 
limit. This is the way how the mean field approach accounts for the exclusion interaction 
of the original process. This, however, not the only way to realize hindrance of the flow 
by the occupancy of the target site. Remaining at the factorized character of the mean 
field current, one could use instead of 1 — p an arbitrary function C(p) of the density of 
the target site for which the following properties are required. First, ({0) = 1, which 
means that, if the target site is empty, there is no hindrance for the current. Second, 
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C(p) is continuous and monotonically decreasing with p and finally ({1) = 0, which is 
responsible for exclusion. The dynamics of this general exclusion model is defined by 
the equations: 

^ = {Pi-ipi-i + qipi+i)C{pi) - [PtC{pt+i) + qi-iC{pi-i)]pi- (36) 

Our aim with the generalization of the original model is to point out that the concrete 
form of C(p) is irrelevant regarding the dynamics of the system in the sense that it 
influences only the random prefactor in the expression of the current through a barrier. 

To see this, the calculations of the previous sections can be carried out with slight 
modifications. With the variable Ui = Pi/C{pi) we obtain in the steady state: 

Vi+i = TiVi - JgrMC(Pi)C(Pi+i)]~^- (37) 
This leads to the same formula as given in Eq. (fT3l) . however, with 

L-l 

Al = e-^- 5^g7^[C(p.)C(p.+i)]-^e^^-. (38) 
i=o 

To obtain an upper bound on we can use inequality f|T5l) which still holds. First, 
let us consider the factors [C(Pj)] ""^ where pj > 1/2. For these factors we can write: 
[({pj)]~^ = Vj/Pj < 2?/j < 2yQe~^K For the factors with pj < 1/2 we can use the 
mo notonicity of C(p) to obtain an upper bound: [C(Pi)]~^ < [C(l/2)]^^ = const. We have 
thus [C(Pi)]~^ ^ max{2|/oe~'^-'', [C(l/2)]~^}. Using these inequalities an upper bound on 
Ax is obtained which contains an 0(1) contribution in the sum for sites where the 
magnitude of the potential is very close to Ul- The backward iteration can be done in 
an analogous way and similar conclusions for A_2,' can be drawn. The final conclusion 
is that the current can be written in the form given in Eq. (1251) in case of an open 
barrier an in the form given in Eq. fl32l) in case of an infinite system and in both cases 
the concrete form of C(p) influences only the prefactors in front of the exponentials. 
Moreover, if relations in Eq. (l26l) are fulfilled, the linear contributions A° and A° ^, are 
independent of the form of C(p) and are thus given by the expressions obtained in the 
previous section. 

5. Phenomenological theory of the dynamics 

The description of the non-stationary state of the system is based on that segments 
of characteristic length ^ can be regarded as quasi-stationary and as time elapses the 
characteristic length scale ^ increases. Thus, the steady state properties of a finite 
system has to be reviewed first. 

5.1. Driven phase 

Let us assume that the system is driven to the right on average, i.e. Af/ < 0. In a 
finite but large system of size iV, there are 0{N) barriers and the stationary current 
J{N) is roughly equal to the smallest one among the currents of barriers J, considered 
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in the previous section. We have obtained there that the current through a barrier 
is Ji = CjC"^'/^ where Ui is the height of the potential barrier and Ci is an 0(1) 
random factor which depends on the shape of the barrier (and that of the environment 
in close vicinity of the barrier). Although we have assumed there that the barrier is 
well separated from other barriers, which does not hold in a disordered system, the 
neighboring barriers are expected to influence only the random prefactor Cj. The 
relevant factor in Jj is e"*^'/^, the inverse of which is roughly the square root of the 
waiting time Tj ~ e^^ of a single random walker at that barrier. The distribution of the 
random variable Xj is known to have an algebraic tail [22] 

P>(r)~Ar~^, (39) 

where the control parameter fi is the positive root of the equation 



e^^^ = 1. (40) 

distribution 



— 1/2 

It follows then that the current through barriers Ji ^ t- has the asymptotic 



P<(J)^yl'j2M (41) 

for J — 7- 0. The current in a finite system then follows a Frechet distribution and has 
the typical value vanishing with as [T5| IT6] 

Jtyp(Ar) ~ Ar-V2^ ^42) 

In the steady state, a phase separation can be observed: almost all mass accumulates 
behind the highest barrier and forms a high density phase of macroscopic size, where the 
density is close to one; in the rest of the system the density is close to zero. Within these 
phases, the density profile is not completely flat but it contains peaks at those barriers 
whose height is greater than the half of the highest potential barrier. The number of 
these peaks is 0{N^^'^) where the exponent 1/2 is universal in the driven phase and is 
related to the half-filling of the highest barrier [T5] . 

Let us assume now that the system is started from a state with random local 
densities and with an average density 1/2. After time t has elapsed, the characteristic 
length scale is ^{t) and the typical size of high and low density segments is ^(t)/2. The 
rate of growth of these domains is proportional to the typical current Jtyp{t) at that 
time scale, that means we can write 

§-Jum-r'^''- (43) 

Integrating this differential equation yields 

^(t) ^ t2A./{l+2M) (44) 

and 

JtyAt)-t-'^^'^'^^- (45) 

So, the typical current decays algebraically just as in the pure model but with a non- 
universal decay exponent /3 = 1/(1 + 2/i). The growth of the length scale ^{t) follows 
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also a power-law. It is, however, more convenient to measure the average distance l{t) 
between adjacent peaks of the density profile in numerical simulations rather than C,{t). 
This quantity grows as 

~ ~ (46) 

again with a non-universal coarsening exponent (5 = /i/(l + 2yu). 

So far we have tacitly assumed that the barriers comprise many sites such that it 
is reasonable to speak of half-filling of barriers. This is, however, not always true when 
the distribution of forward hop rates p is not bounded away from zero. In that case the 
barriers may typically consist of single links through which the rate of forward hopping 
is vary small, and as a consequence, the above calculations have to be modified. Let us 
assume that the distribution of p has the asymptotic form P<{p) — const ■p'^ for p — 0. 
Comparing this to Eq. ( HTj) . it is clear that whenever v < 2ix the local currents are 
controlled almost always by barriers consisting of single links and Eq. ( H2|) changes to 

Jtyp(iV) ~ iV-V- (z. < 2/x). (47) 

This anomalous scaling of the stationary current has been revealed in Ref. |I5j. Here 
we go further and derive how the dynamics are modified if z/ < 2/i. Using Eq. (HTj), 
the evolution equation ^ ~ JtypHif)] for the typical length of quasistationary segments 
results in 

e(t)~r/('+"\ (48) 

This relation together with Eq. (147|) yields the following time dependence of the typical 
current: 

Jtyp(t) ~ (z/<2^). (49) 

In a segment of size ^, where the quasistationary current is J(^) ~ ■C^^^'^, mass 
accumulates at those extended barriers where the waiting time ~ e^' is greater than 
1/J(^). Making use of the distribution of waiting times in Eq. ( 139|) we obtain that the 
number n(^) of such barriers in the segment scales as 

nio-iPAi/jm-e-''". (50) 

Thus, the length scale l{t) grows with time as 

m ~ i/n{i) ~ e'" - t^'^'^''^ {y < 2/x). (51) 

5.2. Zero average force 



If the average force is zero, i.e. AU = 0, then the extension of the largest barrier is 0{N) 
and the above theory breaks down. At this point we have only scaling considerations 
at our disposal [15j. The height of the largest barrier is 0{\/N), therefore the typical 
stationary current in a finite system is expected to scale with N as 



InJtyp(iV) ~ ViV. (52) 
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Plugging this relation into the r.h.s. of Eq. f l43p yields 

e(t)~[ln(t/lnt)f (53) 

and 

Jtyp{t) t-Hnt. (54) 



In the steady state of a finite system almost all mass is concentrated in 0{^/N) basins 
where the density is close to one and the extension of the largest basin is 0{N) |15j . 
Thus the number of jumps in the density profile, where the density crosses over from high 
(p ~ 1) to low (p ~ 0) density is 0{\/N). Defining l{t) as the average distance between 
adjacent jumps in the profile at time t and assuming quasi-stationarity in segments of 
characteristic size ^(t), we obtain that it increases with time as 

VW)~ln(t/lnt). (55) 

In case of zero average force, we have formally /i = from Eq. (j40|) . The formulae 
for the dynamical quantities obtained here are consistent (apart from logarithmic 
factors) with those valid in the driven phase taken in the limit p — )■ 0. 

6. Numerical analysis 

The stationary properties of the disordered model predicted by the phenomenological 
theory has been compared with results of Monte Carlo simulations and a good agreement 
has been found [15]. The dynamical behavior of the current and the length scale has not 
been directly checked. The reason for this is that even for each random sample many 
runs have to be performed with different stochastic histories in order to calculate to local 
currents or the density profile. Instead of this, the time dependence of the displacement 
of a tagged particle has been measured in case of zero average force [9] and in the driven 
phase [15]. By solving the evolution equations ([1]) in the mean field treatment, the local 
densities and currents are directly at our disposal and the dynamical behavior of Jtyp(^) 
and l{t) can be conveniently checked. 

In the numerical calculations, we have considered two types of distributions for the 
hop rates. A discrete one, where pi + g^+i = 1 and the probability density of p is 

f{p) = c6[r/{l + r) - + (1 - c)5[l/(l + r) - p], (56) 

where < c < 1/2 and < r < 1 are constants, and a continuous one with probability 
densities 

f{p) = 1/s if < p < s and f{p) = otherwise; 

g{q) = 1 if < g < 1 and g{q) = otherwise. 

The control parameter fi is given by 

_ ln(c-i - 1) 
^ ~ ln(l/r) 



(57) 
(58) 
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in the former case and implicitly by 

s = (l-^2)-i/M (59) 

in the latter case. In the case of the continuous randomness, the anomalous scaling given 
in Eqs. ( 149|) and ( IST]) sets in with i/ = lif/i>l/2, while for the discrete randomness 
the scaling is never anomalous. 

We have generated random samples of size L = 10^ — 10^ and starting from a 
disordered initial state where the local densities are independent random variables with 
a homogeneous distribution in the range [0, 1], the evolution equations in Eq. ([T]) have 
been numerically integrated by the 4th order Runge-Kutta method [23] up to time 2^^. 
For the times where measurements were carried out, the coarsening length scale was 
much less than the size of the system so that the system can be practically regarded as 
infinite. We have calculated the time- dependence of the finite-L estimate of the typical 
current given in Eq. (jS]) and the time-dependence of the length l{t) = L/n(t), where 
n{t) is the number of points where the density profile crosses the line p = 1/2. These 
calculations have been repeated for 10^ independent random samples and the averages 
of the above quantities have been calculated. Having the measured data Jtyp(^n) and 
l{tn), we have calculated effective exponents from neighboring data points at time tn 
and tn+i- 

^ ln[Jt.p(Wl)/Jtyp(tn)] ( ) ^ HKtn,.) /Ktn)] _ (go) 

In addition to this, we have also investigated the distribution of local currents. 

We start the presentation of numerical results with the driven phase, where /i > 0. 
The distribution of local currents at different times can be seen in Fig. [1] for the binary 
randomness and in Fig. [2] for the uniform one. 




Figure 1. Left: Numerically calculated distribution of the logarithm of local currents 
at different times for the binary randomness with parameters r = 1/4, c = 1/3. The 
control parameter is — 1/2 and (3 = 1/2. Right: Scaling plot of the same data. 



As can be seen in the figures, an adequate data collapsing can be achieved using 
the scaling variable Jt^ where /3 is the exponent predicted by the theory in Eq. ( H5|l . 
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Figure 2. Left: Numerically calculated distribution of the logarithm of local currents 
at different times for the uniform randomness with s = 1.5. The control parameter is 
/i w 0.37607 and /3 w 0.570734. Right: Scaling plot of the same data. 



The time-dependence of the typical current has been calculated in several points of the 
driven phase and the corresponding effective exponents f3esit) are plotted against time 
in Fig. [3l The obtained data are again in satisfactory agreement with the predictions 
of phenomenological theory. 




Figure 3. Effective exponent /3cs{t) plotted against time in different points of the 
driven phase. Left: Data are obtained with the binary randomness with parameters 
c = 0.2, r = 0.5 and c = 1/3, r — 0.5,0.25,0.125 where the control parameter is 
/i = 2, 1, 1/2, 1/3, respectively. Right: Data are obtained with the uniform randomness 
with s = 3,2,1.5,1.29 where /i « 0.7488,0.5722,0.3761,0.25, respectively The 
horizontal lines indicate the value of /3 predicted the phenomenological theory. 



The length l{t) has been measured at the same points of the driven phase, 
as well. The corresponding effective exponents are compared to the predictions of 
phenomenological theory in Fig. HI As can be seen, the finite time corrections are more 
considerable than those of (3, nevertheless the asymptotic behavior is still compatible 
with the theory. 

Next we turn to present numerical results obtained for zero average force (/x = 0). 
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Figure 4. Effective exponent Scs{t) plotted against time in different points of tlie 
driven phase. Left: Data are obtained witfi binary randomness with parameters 
c = 0.2, r = 0.5 and c — 1/3, r = 0.5,0.25,0.125 where the control parameter is 
/i = 2, 1, 1/2, 1/3, respectively. Right: Data are obtained with the uniform randomness 
with s = 3,2,1.5,1.29 where n « 0.7488,0.5722,0.3761,0.25, respectively. The 
horizontal lines indicate the value of S predicted the phenomenological theory. 



The distributions of local currents for different times are shown in Fig. |5] and |6l The 
striking difference compared to the driven phase is that, although the exponent /3 is still 
finite (/3 = 1) for = 0, the distributions are broadening with increasing time. A rough 
scaling collapse can be achieved in terms of the scaling variable ln(J)/ln(t). Earlier 
results on the distribution of the stationary current in finite systems of size L showed 
an approximate scaling collapse for the scaling variable ln( J)/L^/^ \15\. Taking into 
account the relation between time and length scale in Eq. (|53|) . this is consistent with 
our present results on the dynamical scaling. As can be seen in Fig. [TJ the effective 




Figure 5. Left: Numerically calculated distribution of the logarithm of local currents 
at different times for the binary randomness with parameters ^ = 1/4, c=l/2. At 
this point the average force is zero and fi — 0. Right: Scaling plot of the same data. 



exponent /3esit) overshoots the expected asymptotical value 1 by a few percent. This 
is in accordance with that the scaling collapse of distributions is not perfect and the 
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Figure 6. Left: Numerically calculated distribution of the logarithm of local currents 
at different times for uniform randomness with parameters s = 1. At this point the 
average force is zero and fj. — Q. Right: Scaling plot of the same data. 



shape of distributions is still slightly changing at the numerically available time scales. 
We have also measured the variance of local currents which enhances the contribution 
of large local currents compared to the typical value defined in Eq. ([3]). These lie just 
in the still deforming and thus poorly scaling part of distributions. The corresponding 
effective exponents are approaching the theoretical value from below, see Fig. [71 In the 




Figure 7. Left: The effective exponent /3off(i) (open symbols) plotted against time 
for zero average force. Filled symbols show the effective exponents calculated from the 
variance of the current. The solid line corresponds to the form in Eq. ([Ml) . Right: 
Time-dependence of the length scale l{t). The effective exponents are shown in the 
inset, where the solid lines corresponds to the form in Eq. ([55t . 



right panel of Fig. [TJ the measured length scale is plotted against time. As can be seen, 
the data are in good agreement with the law given in Eq. (!55|) . 
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7. Discussion 



The partially asymmetric simple exclusion process with random-force disorder has 
been investigated in earlier studies exclusively by Monte Carlo simulations and by a 
phenomenological random barrier description. In this work, a mean field approximation 
has been applied to this model and the main focus was on the non-stationary phenomena. 
The mean field approximation leads to a system of deterministic, nonlinear differential 
equations. This model is less complex than the original stochastic process but still not 
tractable analytically. Nevertheless, it is appropriate for applying a phenomenological 
random barrier theory to it. According to our analytical mean field calculations, the 
key issue of the theory, namely the current through a barrier has the same behavior as 
that has been assumed intuitively for the original model in earlier works. This leads 
to the same large-scale stationary and non- stationary behavior as has been conjectured 
for the original model. We have investigated the mean field model numerically, which 
is considerably faster than performing Monte Carlo simulations and have found that 
the dynamics is satisfactorily described by the phenomenological theory. Since a good 
agreement between the phenomenology and results of Monte Carlo simulations carried 
out on the original model has been found in earlier works, we conjecture that the mean 
field model belongs to the same universahty class as the original stochastic process does. 
This means that the static and dynamical exponents in the driven phase, as well as the 
scaling relations for the case of zero average force are identical. Our results show that 
in the presence of disorder the local correlations are unimportant concerning the large 
scale behavior of the system. This conclusion can be instructive for the investigation 
of other transport processes with random-force disorder, where the simple mean field 
approximation may give the correct large scale behavior. 



Appendix A. 

Let us denote the stationary local density in the pure model by p and the deviation 
from the stationary density by ei(t), i.e. pi(t) = ei(t) + p. The spatially continuous limit 
of the evolution equations Eq. ([T]) reads as 

de ^ d^e dt ^ de ^ ^ 

dt dx"^ dx dx^ 
where the constants D, v and A are given in terms of the jump rates as D = [p + q)/2, 
V = [1 — 2p){p — q) and A = 2{p — q). The local current in the continuum limit takes 
the form 

J{x,t) = J^ + ve-D—-Xe\ (A.2) 

where Joo = {p — '?)p(l ~ p) is the current in the steady state. First, let us consider the 
simplest case p = q, which describes the symmetric simple exclusion process. In this 
case, V = X = and Eq. (lA.lll reduces to the diffusion equation. Consider a random 
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initial density profile of the form 

oo 

e(a;,0) = ^ sj{x-n), (A.3) 

n=— oo 

where 6{x) is the Dirac delta distribution and the s„ are independent binary random 
variables with the probability density /(s) = ^6{s — 1) + ^S{s + 1). The solution of Eq. 
(lA.ip is then 



1 v-^ —(x — n — vtV . 



The mean value of local quantities such as the deviation e(t) = limi;,^oo 2Z+1 '^n=-L ^(^' ^) 
at time t can be calculated alternatively from e{t) = J e(0, t)f{s)ds since the Sn are iden- 
tically distributed for all n. For the mean deviation we obtain the obvious result e{t) = 0, 
since e(0) = and the total mass / p{x, t)dx is conserved by Eq. (lA.ip . The fluctuations 
of e{x, t) are characterized by the variance, the square of which can be easily calculated: 

i^M = -J— Y exp ^ ^ -—- / exp ^dx = {MDty^/\ (A.5) 
^ ' A-nDt ^ ^ 2Dt A-nDt J ^ ^ 2Dt ^ ' ' ^ ' 

n=-oo ^ 

where the sum has been approximated by an integral. Thus, the typical deviation from 
the stationary density measured at a randomly chosen site at time t is in the order of 
^-1/4 rjj^g typical local current can be estimated in a similar way. The mean value ^ 
is zero while the square of the variance is 



O \ 2 ^00 9 2 

oe \ 1 n —n 



- 



E 



dx ) A-nDt ^ ADH'^ 2Dt 

The typical current measured at a given bond at time t is thus 0{t~^^'^). 

In the case p q, which corresponds to the asymmetric simple exclusion process, 
the Galilean transformation x ^ x + vt cancels the second term on the r.h.s. of Eq. 
f lA.ip and one obtains the noiseless Burgers equation 



de ^ d'^e \ ^7 
dt dx"^ dx 

This can be exactly solved by the Cole-Hopf transformation (see e.g. [24J) w{x,t) = 
exp 2^ e{x' ,t)dx' , which maps the Burgers equation to the diffusion equation. The 
solution is obtained from 

2D d\nw , , ^, 



with 



w{x,t) = I ^y^^^^P [-^-^J^ + ^ I e{x",0)dx"]dx'. (A.9) 
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Using the random initial condition given in Eq. (lA.3p . the deviation of the density at 
X = is given as 

Dt J_^exp[-f^^ + ^r{x)\dx 

where r(x) = e{x' ,0)dx' is a piecewise constant function with unit jumps at integers 
whereas for non-integers it is given by r{x) = — ^ + X^nio'^™' where [x] denotes the 
integer part of x. Thus r(x) can be regarded as a random walk which makes jumps at 
integer "times" x. It is easy to see that the mean value of the deviation e{t) is zero since 
exp 2^r(x) is an even function of x due to Prob[r(x)] = Prob[r(— x)]. The typical value 
of the deviation at time t can be obtained as follows. First notice that the r.h.s. of Eq. 
( lA.lOp can be regarded as the expected value of x which has the (unnormalized) weight 



function expF(a:) with 

Fix) = + —r(x). (A.ll) 

^ ' ADt 2D ^ ' ^ ' 

The dominant contribution to this expected value comes from the interval where F{x) is 

maximal, since otherwise the weight exp F{x) is negligible. So we have the approximate 

relation e(0,t) ~ -^x^aax, where Xmax is the location of the maximum of F{x). (In 

case there are many maxima, Xmax is their mean value.) The distribution of Xmax is 

symmetric around zero and the dependence of its magnitude on time can be obtained by 

taking into account that the variance of r(x) which characterizes its typical magnitude 



is yr'^{x) = \/[x] + 1/4, i.e. proportional to ^/x for large x. Replacing r{x) in Eq. 
( lA.llI) by ^/x, we obtain a non-random function the maximum of which is at {2Xt)'^^^. 
Thus the width of the distribution of x^ax is in the order of t^^^ and we obtain finally 
that the typical deviation from the stationary density measured at a given site scales 
with time as etyp(t) ~ t~^^^. By a similar calculation one can show that the expected 
value of ^ is of the order of t^^^. Now we can turn to the analysis of the fluctuations 
of the current. If p 7^ 1/2 then f 7^ and the fluctuations are dominated by the term 
ve, see Eq. (]A.2|) . Thus the magnitude of the typical local current (relative to the 
stationary current) scales as (J — Joo)typ(^) ~ t^^^^. If, however, p = 1/2, the above 
term is zero and the fluctuations are determined by the other two terms leading to 

(J- Joo)typ(t)~t-'/'. 
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